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Measurement of the high energy y-rays from heavy ion reactions using Cerenkov detector” 
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The energetic bremsstrahlung photons up to 100 MeV produced in heavy ion collisions can be used as a 
sensitive probe for short-range correlation in atomic nuclei. The energy of the y-rays can be measured by col- 
lecting the Cerenkov light in the medium induced by the fast electrons generated in the Compton scattering or 
electromagnetic shower of the incident y ray. Two types of detectors based on pure water and lead glass as sen- 
sitive materials were designed for this purpose. The y response and optical photon propagation in the detectors 
were simulated based on electromagnetic and optical processes in Geant4. The inherent energy resolutions of 


0.022(4) + 0.51(2)/E}/? for water and 0.0026(3) + 0.446(3) / E+’? for lead glass were obtained. The geom- 
etry sizes of the lead glass and water were optimized to 30 cmx 30 cmx 30 cm and 60 cmx 60 cmx 120 cm, 
respectively, to detect high-energy y-rays at 160 MeV. The Hough transform method was applied to reconstruct 
the direction of the incident y-rays, providing the ability to experimentally distinguish the high-energy y-rays 
produced in the reactions on the target from random background cosmic-ray muons. 
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I. INTRODUCTION 


Bremsstrahlung high-energy photons produced in heavy- 
ion reactions have attracted increasing interest because of 
their relevance to the nuclear equation of state (nEOS) and 
their short-range correlation in nuclei. In nEOS studies, 
particularly for nuclear matter with large neutron-to-proton 
asymmetry, a variety of isospin probes have been identified 
to constrain Esys(p) (the density-dependent nuclear symme- 
try energy), including the preequilibrium n/p yield ratio [1], 
n/p differential flow [2, 3] and bremsstrahlung high-energy 
photons [4]. Among these probes, bremsstrahlung y-rays cre- 
ated by heavy ion collisions are clearly observable because 
of their rare interactions with the medium after they are pro- 
duced. Very recently, it has been revealed that bremsstrahlung 
high-energy -y-rays carry information on the high-momentum 
tail (HMT) of nucleons, giving rise to the short-range corre- 
lation effect in nuclei [5-8]. On the other hand, however, the 
experimental data in this direction is quite scarce. 

Recently, the full y energy spectrum up to 80 MeV was 
measured in reactions ®°Kr+!?4Sn at 25 MeV/u with a 15- 
unit CsI(T1) hodoscope mounted on a compact spectrometer 
for heavy-ion experiments (CSHINE) [9-13]. It has been 
demonstrated that the y energy spectrum above 20 MeV is 
reproduced fairly well by transport model simulations that 
incorporate ~y production from incoherent np scattering with 
an approximate 15% HMT ratio [14]. However, CsI(TI) is a 
slow detector, and the microsecond response time of CsI(T1) 
crystals makes it difficult to reconstruct the total energy from 
multiple firing units. Therefore, we are motivated to develop 
a fast and relatively cheap detector containing a sufficiently 
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large-volume-sensitive material to detect high-energy rays 
in heavy-ion reactions. The Cerenkov radiation [15] detec- 
tor is a favorable option because of its fast response time in 
the order of tens of nanoseconds and its ability to infer the 
incident direction information of the initial y-rays, the latter 
of which can be used to suppress the cosmic-ray muon back- 
ground from random directions. 


In this paper, we report the design of a Cerenkov y 
calorimeter using water and lead glass as sensitive media. 
Based on Geant4 packages, the geometric size of the detectors 
was optimized. The energy resolution was obtained by track- 
ing each Čerenkov photon before it arrived at the photomulti- 
plier tube (PMT) for which the quantum response was mod- 
eled. The incident direction reconstruction was implemented 
using the Hough transform method. The remainder of this 
paper is organized as follows. Section II describes the simu- 
lations framework of the calorimeter. Section II presents the 
optimization of the detector size and the reconstruction of y 
direction. Finally, Sect. IV concludes the paper. 


Lead glass 


Fig. 1. (Color online) Detector configuration with two sensitive vol- 
umes, lead glass (left) and pure water (right), respectively. 
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Il. SIMULATION SETUP 


In this study, Geant4 (version 4.10.05)[16] packages were 
used for Monte Carlo simulation and optimization of the de- 
tector. “QBBC” and “G4OpticalPhysics” are applied as the 
physical process list to describe the electromagnetic (EM) 
showers of y rays in materials, and to model the generation 
and transport of Cerenkov photons. For each event in the sim- 
ulations, incident y-rays hit the front of the detector. Then, 
Cerenkov photons are generated if fast electrons are produced 
by Compton scattering or e~e* generation. Each Cerenkov 
photon is tracked to its termination, either to be absorbed dur- 
ing propagation or to reach the surface of the PMTs, where 
the waveform pulse of the given parameters is generated with 
a certain quantum efficiency. The waveforms were recorded 
at intervals of 2 ns for digitization. The final data correspond- 
ing to each incident y ray are saved as a matrix of N; x N2 
dimensions, where NV; represents the number of fired PMTs 
and Nz represents the number of sampling points for the cor- 
responding waveform. 


A. Detector geometry 


The detector structure and locations of the PMTs are shown 
in Fig. | for lead glass and water as sensitive materials, de- 
fined as “G4_GLASS_LEAD’ (left) and “G4_WATER’ (right), 
respectively. The water tank size was 60 cm x60 cm x 120 cm, 
and the size of the lead glass was 30 cmx30 cmx30 cm. The 
PMTs were arranged in an 8 x 8 array with a water configura- 
tion. In the lead glass configuration, the PMTs were arranged 
in 4 x 4 arrays on the four sides of the detecting tank. The 
diameter of each PMT was 51 mm, and the distance between 
each neighboring PMT pair was 70 mm, both vertically and 
horizontally. 


B. Optical process 


Upon invoking the Cerenkov mechanism in Geant4, the en- 
ergy and number of Cerenkov photons were sampled in each 
G4step according to[17]. 
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where 6, is the Cerenkov angle, \ is the wavelength of the 
Cerenkov photon. The initial position of Cerenkov photons is 
uniformly distributed in every G4step, the emission angle is 
calculated according to the refractive index of materials and 
the speed of the charged particle, the outgoing azimuth is uni- 
formly distributed within the range of 27, we set a maximum 
of 100 photons emitted in each step to ensure the detailed 
sampling. In the process of photon transport, the transmis- 
sion characteristics of photons in a material and their behavior 
at the boundary between the two materials must be defined. 
In this simulation, we defined the scattering and absorption 
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lengths between the Čerenkov photons and water molecules 
by referring to the test data of IceCube [18, 19]. Owing to the 
lack of optical parameters for lead glass, we conservatively 
defined an attenuation efficiency of 70% for a 10 cm prop- 
agation. For the boundary characteristics [20, 21], we used 
UNIFIED model [22, 23] in Geant4 and selected “dielectric- 
dielectric” option to describe the interface between the mate- 
rial and PMTs. In this model, Geant4 determines the photon 
boundary behavior according to the Fresnel formula and the 
refractive index on both sides. For the remaining boundaries, 
we use the dielectric- LUT model [24] and select a polished 
Teflon -LUT boundary. Thus, Geant4 determines the reflec- 
tion, refraction, and absorption of photons based on the built- 
in parameters. 


C. PMT Response 


In the full case, photons are converted into photoelectrons 
with a certain quantum efficiency after hitting the PMT, and 
a pulse is formed after multiplication. A pulse formed by a 
single photon is described in [25] 
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where ti = thit + ttrans and thit represent the time at which 
a photon hits the PMT, tirans = 29 ns represents the elec- 
tron transit time of the PMT, o = 1.2 ns represents the transit 
time spread. The final waveform is generated by superimpos- 
ing all single-photon waveforms when multiple photons are 
converted into photoelectrons, as shown in Fig. 2 (a). Based 
on the incidence of y rays, the waveform of each PMT within 
240 ns was recorded as the final data. Figure 2(b) shows the 
distribution of the time at which the optical photons reach the 
PMTs in the lead-glass configuration, which was extracted 
by linearly fitting the rising edge of the waveform [26]. This 
illustrates that most photons reach the surface of PMTs be- 
tween 26 ns and 31 ns after y emission, which means that 
we can distinguish between direct and scattered photons ac- 
cording to the distribution in the direction reconstruction (see 
Sect. IID). 
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Fig. 2. (Color online) (a) a typical waveform for PMT, (b) the dis- 
tribution of the time when optical photons reach PMTs in lead glass 
detector. 
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Ill. RESULT AND DISCUSSION 


A. Influence of detector size on energy resolution 


We used the photoelectron peak number (npe) and the en- 
ergy resolution, defined by 6g, = dbn,. = Onpe/ (Npe), to 
optimize the detector design. In the simulation, such high- 
energy y rays hit the center of the front surface of the de- 
tector perpendicularly. The shower electrons and positrons, 
if produced with a velocity exceeding the speed of light in 
the medium, will generate Cerenkov light propagating to the 
PMT, where the photoelectrons are generated. Owing to the 
statistical fluctuations, the number of photoelectrons varies. 
Figure 3 (a) presents the distribution of the number of photo- 
electrons for 50 MeV incident y-ray in the lead-glass detector 
as an example. In the following analysis, the photoelectron 
peak number (npe) was considered as the average number of 
photoelectrons. Figure 3 (b) presents the distribution of (npe) 
as a function of incident y energy E- for the two configura- 
tions at their own optimized volumes (see below). For a de- 
tector of a given size, (Npe) has a linear dependence on Æ}. 
Thus, the y-ray energy can be measured using the number of 
photoelectrons. 
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Fig. 3. (Color online) (a)the spectrum of the photoelectron yield 
for 50 MeV y rays in lead glass detector,(b) Linear response of the 
calorimeter. 


We then optimized the detector size at a given maximum y 
energy of 160 MeV, which covered the range of interest for 
E, in heavy-ion reactions at Fermi energies. For each event, 
the distribution of photoelectron number was analyzed to ob- 
tain (npe) and its standard deviation (on,,). If the detector 
medium is too small, much of the y-ray energy will leak out- 
side the sensitive volume. However, if the detector medium 
is too large, Cerenkov photons are scattered many times and 
gradually absorbed, leading to a reduction in the number of 
photoelectron collected by PMTs. These two factors collec- 
tively determine the energy resolution. Figure 4(a) and (c) 
illustrate the energy resolutions ôg and (npe) as functions of 
the horizontal and longitudinal lengths of the lead-glass con- 
figuration. Clearly, as the horizontal and longitudinal lengths 
increased, 6 E, Teaches its lowest point at 30 cm, whereas 
(Npe) decreases when the horizontal or longitudinal lengths 
exceed 30 cm. Thus, 30 cmx 30 cmx 30 cm is the optimal 
size of the lead glass. Figure 4(b)(d) show the same quantities 
for the pure water configuration. As the longitudinal length 
increases, the energy resolution ôg, decreases gradually and 
converges to 6%, and (npe) first increases and then decreases 
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because Cerenkov light attenuation becomes the main factor 
after the longitudinal length exceeds 80 cm. As the horizontal 
length increases, ôg, reaches its minimum at 60 cm, beyond 
which (npe) starts to decrease. Thus, 60 cm x 60 cm x 120 
cm was the optimal size for water. 
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Fig. 4. (Color online) process of size optimization.(a)relationship 
between energy resolution and detector size for lead glass, 
(b)relationship between energy resolution and detector size for wa- 
ter, (c) relationship between (npe) and detector size for lead glass, 
(d)relationship between (npe) and detector size for water. 
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Fig. 5. (Color online) The initial(a) and reconstructed(b) y energy 
spectra in lead glass configuration. 


Given the good linear response of the water and lead glass 
Cerenkov calorimeter to the y energy, as shown in Fig. 3 
(b), one can reconstruct the y energy from the signal height 
equivalent to the number of photoelectrons. To test this abil- 
ity, we simulated the detector response for 10° y events with 
an initial energy Finitial in an exponential distribution. The 
slope of the input exponential distribution is set to —0.05, 
as shown in Fig. 5 (a). The reconstructed energy (Erec) 
is plotted in panel (b) with the slope parameter fitted at 
—0.049. It is shown that the Cerenkov calorimeter of lead 
glass measures high-energy y in the range from 5MeV to 
160MeV. Figure 6 shows the resolution at various incident 
energies for the lead glass and water configurations. Inher- 


ent resolutions of 0.022(4) + 0.51(2)/Ey/” for water and 


0.0026(3) + 0.446(3) / Ey’? for lead glass were obtained by 
fitting the simulated data points. At high energies(above 100 
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Fig. 6. (Color online) Resolution prediction of the calorimeter of 
water and lead glass, respectively. 


MeV), the resolutions were saturated at approximately 7.3% 
and 4.7%. 


B. Direction reconstruction 


It is well known that a definite angle exists between the 
Cerenkov photons and charged particles [17], which is the 
basis for direction reconstruction. In fact, ~y shower also par- 
tially retains this feature. Figure 7 shows the angle distribu- 
tion between the Cerenkov radiation direction and the initial 
direction of electrons (y rays) in water (lead glass), which 
was obtained by the Geant4 simulation, where the energies 
of electrons and y were sampled evenly from 5 to 160 MeV 
in the simulation. The refractive index of lead glass and 
water are 1.7 and 1.3, so the cosine of their Cerenkov an- 
gle are cos@, œ~ 0.58 and 0.77, respectively. According to 
Fig. 7, although the e*e~ pair production and Compton ef- 
fect may cause scattering, the emission angle distribution of 
the Cerenkov photons produced by the EM shower is still re- 
lated to the initial direction of y rays. This suggests that the 
direction of y rays can be reconstructed by referring to the 
electron direction reconstruction method used in large experi- 
ments, such as the Super-Kamiokande and Sudbury Neutrino 
Observatory(SNO) [28, 29]. It was found in our work that 
the Cerenkov photons experience scattering many times be- 
fore reaching PMTs in water because of the overlength of the 
medium, heavily smearing the initial direction information; 
therefore, we only reconstructed the y ray direction in the 
lead glass configuration. 


1. Vertex reconstruction 


To reconstruct the direction of the electrons, it is usually 
assumed that the electrons emit Cerenkov light from a fixed 
point. According to the angle distribution in Fig. 7, it can 
be assumed that the 7 rays emit Cerenkov light from a fixed 
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Fig. 7. (Color online) Čerenkov photon direction distribution for 
electron incidence in water and y ray incidence in lead glass, re- 
spectively. 


point with a specific Čerenkov angle; hence, the time taken 
for photons to reach the PMT can be expressed as [30, 31] 


— | Xoi T Xna | 


thit = + to (3) 
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where to represents the moment when the Čerenkov light is 
generated, thi represents the moment when the photon hits 
PMTs, v is the velocity of light in lead glass, Xpmt and 


=> 


Xvtx are the coordinate of the PMT and the vertex respec- 
tively. In our analysis, the optimal estimation of the ver- 
tex coordinates is obtained by minimizing the x? of fitting 
the time distribution with formula (3), in which the tọ and 
X tx are fitting parameters. For each y event, the timing of 
each PMT was extracted by linear fitting to the rising edge of 
the waveform, where the crossing point of the linear fitting 
and the zero baseline was taken as the timing signal of the 
PMTs [26]. The spatial coordinates of each firing PMT are 
used as X oint Because the reflector layer is set in the simula- 
tion, some Čerenkov photons are reflected before hitting the 
PMTs and the timing signals deviate from (3). So according 
to Fig. 2(b), we only selected the PMTs with the hit time be- 
ing less than 1.5 ns before the peak and 1 ns after the peak of 
the time distribution. We define the coordinates of the center 
of the lead glass as (0 cm, 0 cm, 0 cm). Figure 8 shows the 
x? distribution of the vertex coordinate fitting for a 10 MeV 
y-ray incidence, indicating the optimal vertex coordinate at 
(1.29 cm, —1.11 cm, —6.34 cm). 


2. Hough transform 


The Hough transform [32-34] has been successfully ap- 
plied to identify Cerenkov rings that can map the vector 
space of the vertex-to-PMT direction to that of the electron 
incident direction. An example of this application is Su- 
perKamiokande [35]. Similarly, we can define the vector 
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Fig. 8. (Color online) The x? distribution contour in the coordinate 
space of the vertex fitting. 


from the vertex of the y ray to the firing PMT and the initial 
direction vector of the y ray as V, and V, respectively, where 
0 represents the angle between the two vectors. The probabil- 
ity distribution of 0 is indicated by the red line in Fig. 7. The 
vector space of the incident direction of y-ray was divided by 
100 x 100 according to (cos 6,¢), and the weight of each cell 
can be expressed as 
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Y= f (cos biik), 


1 


= cos Dijk = Vij : Vox (4) 


266 


Vii represents the central unit vector of the cell at row 
i and column j in the vector space of incident direction of 
y ray, Vox represents the unit vector from the vertex point- 
ing to the center of kt? firing PMT, and function f rep- 
resents the Cerenkov angle distribution function of y rays 
(the red line in Fig 7) in lead glass. Figure 9 shows the 
event display of Hough transform for an incident y event 
in the direction (0.068, 0.063, —0.995). Figure 9(a) shows 
the hit PMTS position distribution for this event, where the 
27 color represents the signal amplitude in the corresponding 
277 PMT. Figure 9(b) shows the result of the Hough transform 
27 for the 1% PMT, Fig 9(c) shows the cumulative result of 
279 the Hough transform for all firing PMTs, and the brightest 
230 point in Fig. 9(c) represents the optimal estimate of the inci- 
21 dent gamma direction. For example, the optimal estimates 


267 


268 


269 


270 


27 


272 
273 
274 


27 


a 


vertical number 
f- 


0 
horizontal number 


10 


-0.5 0 


cos 


0.5 


Fig. 9. (Color online) Event display of Hough transformation. (a) 
is the position distribution of firing PMTs, there are four sides to 
place PMTs in lead glass configuration, top, left, right, and back, 
(b) is the result of Hough transform for the marked 1°* PMT on the 
back surface, (c) is the cumulative result of Hough transform for all 
firing PMTs in the time window. The cross indicates the optimized 
direction. 
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Figure 10 shows the distribution of cos A0, where A9 is 
the angle between the reconstructed direction and the initial 
237 direction of Ņ-rays hitting the front of the lead glass uniformly 
288 from the target. The peak of the cosine values is close to 
239 cos AM = 1, indicating that the detector can reconstruct the 
290 direction of the signal in the lead-glass configuration. How- 
21 ever, the cosine distribution broadens considerably because 
292 of the rough assumption that y rays emit Cerenkov light at a 
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Fig. 10. (Color online) The distribution of the angle between initial 
and reconstructed direction of y rays 


fixed point. In fact, according to the red line in Fig. 7, most 
y rays would generate Cerenkov light in a path whose length 
is comparable to the detector size, which contributes to the 
bias in the direction reconstruction. The antisymmetry of the 
locations of the PMTs also causes bias. 


Fig. 11. (Color online) Two-dimensional distribution of cos © and 
(Npe) for y rays(a) from the reaction target and for cosmic ray 
muons (b). 


3. Discrimination between y and cosmic ray muon 


In a real beam experiment, only y rays from the reactions 
on the target are of interest. Because the direction of y rays 
from the reaction target is different from that of the cosmic 
rays, it provides a way to suppress the background. To test the 
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ability to suppress the cosmic-ray background, we generated 
y rays with energies between 5 — 160 MeV at the front of the 
detector and mixed them with uniform u~ emissions from the 
top of the detector. The u~ energy E, (in GeV) and zenith 
angle 0, were sampled using the Gaisser formula [36] 
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Considering that the threshold for ~~ to produce Cerenkov 
radiation in lead glass was 78 MeV, we set the sampling range 
to 80 — 1000 MeV. The physical quantity © denotes the an- 
gle between the reconstructed direction and vector from the 
reaction target to the fitted vertex vector. Figure 11 shows the 
two-dimensional distributions of cosO and (npe) for y rays 
from the target (a) and cosmic rays (b), respectively. A very 
different feature between the reaction y rays and the cosmic- 
ray muon background is evident. The cos © of y rays is con- 
centrated above 0.5, (npe) is relatively evenly distributed be- 
tween 0 and 1000, whereas the cosO of u~ is concentrated 
between 0 and 0.5, and (npe) is concentrated around 200. 
Therefore, the directivity of the Cerenkov light provides new 
dimensional information to distinguish the signal from the 
background. 


IV. CONCLUSION 


In this study, we investigated the feasibility of using a 
Cerenkov calorimeter to detect bremsstrahlung ~y rays from 
heavy-ion reactions at Fermi energies. A full framework was 
established to simulate the response and performance of the 
Cerenkov gamma calorimeter based on Geant4 packages, in- 
cluding y-induced EM shower, Cerenkov photon generation 
and propagation, and the parameterization of PMT waveform. 
The optimal volume, linear response, and energy resolution of 
the detector were obtained using water and lead glass as sen- 
sitive media. The inherent energy resolutions at 0.022(4) + 


0.51(2)/ E1? level for water and 0.0026 (3) +0.446(3)/ E1? 
level for lead glass were predicted. It was demonstrated that 
the initial direction of y rays can be reconstructed using the 
vertex fit and Hough transform method, showing the ability 
to distinguish the bremsstrahlung y rays produced in the re- 
actions from the cosmic-ray muon background in a real beam 
experiment. The detector is built and applied briefly to mea- 
sure the high-energy ~y rays produced in heavy-ion reactions. 
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